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Abstract 

In this paper we study the household-structure SIS epidemic spreading on general 
complex networks. The household structure gives us the way to distinguish inner 
and the outer infection rate. Unlike household-structure models on homogenous 
networks, such as regular and random networks, here we consider heterogeneous 
networks with arbitrary degree distribution p(k). First we introduce the epidemic 
model. Then rate equations under mean field appropriation and computer simula- 
tions are used here to analyze our model. Some unique phenomena only existing in 
divergent network with household structure is found, while we also get some simi- 
lar conclusions that some simple geometrical quantities of networks have important 
impression on infection property of infectous disease. It seems that in our model 
even when local cure rate is greater than inner infection rate in every household, 
disease still can spread on scale-free network. It implies that no disease is spreading 
in every single household, but for the whole network, disease is spreading. Since 
our society network seems like this structure, maybe this conclusion remind us that 
during disease spreading we should pay more attention on network structure than 
local cure condition. 
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1 Introduction 



The spread of disease has been one of the focuses in the field of statisti- 
cal physics for many years. The dynamical behavior of so-called susceptible- 
infected-susceptible (SIS) model and susceptible- infected- removed (SIR) model 
have been widely investigated on regular network and complex networks [1-12]. 
Within the studying, individuals are modeled as sites and possible contacts 
between individuals are linked by edges between the sites. It is easy to see 
that both the properties of disease and topological character of network de- 
termine the dynamics of the spread of disease. Studies have showed that there 
is an epidemic threshold A c on regular networks. If the effective spreading rate 
A>A C , the infection spreads and becomes endemic; otherwise the infection will 
die out. While the threshold disappears on scale- free networks [4]. 

Usually, infectious diseases, such as HIV and computer virus, have the simi- 
lar spreading property. They not only can spread in one household, but also 
can spread from one household to another. To study this spreading charac- 
ter, there have been of considerable interests to epidemic models spreading 
among a community of households [12- 17]. These studies were concerned with 
SIR model, which cannot appear endemic behavior. In 1999, Ball introduced 
the SIS household-structure model [18], in which the population is partitioned 
into m households with N members in each household. A threshold param- 
eter i?* was defined. It is shown that for the household with 2 members, 
if R* < 1 then the epidemic die out; if R* > 1 the epidemic will exist at 
an endemic equilibrium. This model has also been studied on homogeneous 
network by the mean of self-consistent field[19,20]. The similar results have 
been obtained. These previous studies about household-structure epidemic 
model were mainly on regular networks. However, studies have showed that 
a large number of systems, such as Internet, world-wide- web, physical, bio- 
logical, and social networks, exhibit complex topological properties [2 1-23]. In 
particular, small- world properties [24] and scale-free degree distributions [25] 
appear in many real network systems. In this paper, we will analyze the SIS 
household-structure epidemic model on complex networks. The outline is as 
follows: 1) introduction; 2) description of the model; 3) mean-field equations; 
4) steady-state solutions; 5) simulation; 6) summary. 



2 Model 



In complex networks with degree distribution p(k), which is the probability 
that a given site has k connections (links) that connect it with other k sites 
(We say that the given site' degree is k.), there are N individuals that are 
grouped as a household on every site. We assume that these N individuals 
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contact each other fully. A healthy individual may get infected from within 
the household and from outside its household. The parameters A and j3 are the 
infection rates from outside and from within the household respectively. We 
give each site x a number i(i e [0, N]), which means that there are i infected 
individuals in the household at site x. The number of infected individuals at 
a given site x changes according to the following transition rates: 



— > 1 at rate A J2( x , y ) h 

1 — > % + 1 at rate if3 for 1 < % < N — 1 

i + 1 — ► % at rate 7 for < % < N — 1 

In the above expressions (x, |/)means that site x and site y are nearest neigh- 
bors, and we suppose there is a (connection) link between them. Infected 
individuals may infect healthy individuals in their household with rate /3,and 
also can infect healthy individuals in their nearest neighbors with rate A. We 
assume that once a site is infected, infections within the site are much more 
likely than infections from outside, so we can neglect the latter. And also, an 
infected individual in a site can recover with rate 7. We suppose that all the 
individuals in a household have the same external connectivity and do not 
take the birth and death into account. 



3 Mean-field equations 

We now solve the above model with mean-field method. Let u k ^ be the density 
of individuals whose household has % infected individuals and the correspond- 
ing site's degree is k, which means this site has k nearest neighbors. It is worth 
noticing that Y^f=o u k,i — 1- According to the transitions rate described in the 
above section, the evolution equations of w^are written as below[4,6]: 



7u fc) i - \kO k (t)u kfi (1) 

XkQ k (t)u kfi - /3ttfc,i + 7«fc,2 - 7"fc,i ( 2 ) 

(i - l)/3u fcji _i - ifiu k ,i +7«fc,i+i - lUk,i (i e [2, N - 1]) (3) 

(N - l)/3M fc)7V _i - 7M fc)7V (4) 



du kt0 (t) 

dt 
du k ,i(t) 

dt 
du k ,i(t) 

dt 
du k>N (t) 

dt 
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In equations(l)-(4), @k{t) is the probability that a link from a site points to 
another site with at least one infected individual. And the expression of 
is: 



N 

e*W = E#'A)EwjW ( 5 ) 

k' 3=1 

where p(k /k) is the probability that a site with k degrees points to another 
site with k degrees. For uncorrected networks the expression of p(k /k) is [7]: 



Substituting (6) to (5), we get ©&(£) = 0(t) independent of fc: 

= ©(*) = 7^E^')Ei^W (7) 

W k' 3=1 



Now we are going to get steady solutions of Eqs.(l)-(4) 



4 Steady-state solutions 



4.1 For N=2 



The evolution equations are simplified as: 



— |± — = \kQ(t)(l - u k>1 - « fc)2 ) + 7Mfc, 2 - /5«fc,i - 7«fc,i (8) 
— g- t — =/3«m -7Mfc,2 (9) 



In Eq.(8), we have used the equality: w^o + Wfc,i + Wfc,2 = 1- When t — > oo, 
= and = 0, then the steady-state solutions are: 



A 7 A;0 . . 

MM = 7 2 + (7 + ^)AA;0 (10) 
A/%© 

Mfc ' 2 - 7 2 + (7 + /5)AA;0 (U) 
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Substituting (10)and (11) to (7), we get the self-consistent equation of 6: 



Clearly, O = is a solution of Eq.(12), which implies that 1^,0 = 1, = 
and -Ufc i2 = is a steady-state solution of Eqs.(8) and (9). A nonzero steady- 
state solution (That is: u k:i ^ 0, for % > 0) is obtained when 7, (3 and A 
satisfy the following inequality: 



Then we can get the spreading threshold: 

A - ^ {k) (13) 

where (k) = J2kkp(k)i (^ 2 ) = J2kk 2 p(k). In other words, the disease will die 
out when A < A c ; otherwise the disease will pervade the system. Clearly, the 
threshold A c is the function of /3, 7 and ^33: • So the degree distribution of 
networks plays an important role on A c . 

For p(k) = Sk,k c , the network is homogeneous and A c = 7 ^ 2) g a~- We can u ^ 
A c to prevent infection in terms of increasing the recover rate 7 or decreasing 
the site degree k c . 



For p(k) = Ck v (v e (2,3]), the networks are scale-free [19]. When k — > 00, 
f>-0,the 
of A the infei 
SIS model [4] 



-n^r — > 0, the threshold is absent. This fact implies that for any positive value 
of A the infection can pervade the system, which is the same as the standard 



4.2 N>2 



Let -^i = 0(i = 1, 2, • • • N). Suppose U fc =(w M , u fc>2 , • • • , M fc ,v) T and V =(1, 0, • • • , 0) T . 
Considering X^Lo u fc,.j = l,then Eqs.(l)-(4) can be written as: 

SXJ k = - XkQV (14) 



The matrix S is: 
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( -7 - P ~ AA;6 7 - AA;0 -AA;0 
P -2/3-7 7 

2/3 -3/3 - 7 


















-AA;0 





-A&;©\ 





-(7V-l)/3-7 7 
(7V-l)/3 -7 / 



(15) 



Since det(S) = (-7)* - AA;0 £ =i (j - 1)!(-t) jv_j "(-W~ 1 ^ 0, so S" 1 exists. 
Thus: 

U fc = -AfcOS^V (16) 



AT AT 

= n - U fc = (-AJfee)nS _1 V = -AA;0n(S^ 1 V) = -AA;0 ]T n^ 1 

3=1 3=1 

(17) 

where n =(1, 2, • • • , N) and ix, = j 
From (15), we get Sjj 1 : 

s -i = (-7) JV - J '(-^- 1 (j - 1)' (18) 

51 (-7)" " AA;©Ef =1 0' - ^K-^^X"^- 1 
Substituting (18) and (17) to (7), we get the self-consistent equation of 6: 

(-^-A^eEtiO'- 1 ) 1 ^)^-^- 1 



That is: 



e = _J_/v A*'(^)^(-^-'jie \ 



Obviously, = is a solution of Eq.(20). In addition, a non-zero solution 
with 9^0 and Uk,i ^ (i — 1, 2, • • • , N) is allowed if the following inequality 



6 



holds: 




That is: 




> 1 (21) 



j\ > 1 (22) 



From (22), we get the epidemic threshold: 

A - 1 ^ (23) 

where f(N, f3, 7) = X/Jli ^ j', and f(N, /3, 7) is an increasing function 

of iV and 7, but a decreasing function of the recover rate 7. So the epidemic 
threshold is determined by three parameters(iV, f3, 7) and the networks degree 
distribution p(k). We notice that the expression (23) involves multiplication 
of the well-known term jj^r [2,4,6,9], which is closely related to the "average" 
number of secondary infections[7,8]. Not surprising, this result is the same as 
that of the standard SIS model[4]. 

For p(k) = 5k,k c i t ne network is homogeneous. Then A c = ^ -p, we can 
increase the recover rate 7 or decrease the site degree k c and the size of the 
household iV to lift A c to prevent the infectious disease from spreading. For 
large N the threshold is very small. 



For p(k) = Ck v (v G (2,3]), the network is scale-free[21]. When k — > 00, 
jj^j — > 0, then A c = j(]v7J^) (Py ^ 0. So the threshold is absent for scale- 
free network. This implies that for any positive value of A, the infection can 
pervade the system even with high recover rate. 



5 Simulation result 



In above section, we have given the analytical result of the SIS model with 
household structure. We find that for regular network there is an epidemic 
threshold A c ; while for scale-free network the threshold disappears. For com- 
parison, we simulate the model on regular network(see Figl) and on scale free 
network(see Fig2) respectively. For simplicity (without lack of generality), we 
set 7 = 1, N = 4. In Fig.l, we plot the fraction of infected individuals in the 
stationary state, p, for different values of (3 on regular network with k c = 4. 
Obviously, there is a threshold A c for each /3. For f3 = 0.6, A c is 0.026, in agree- 
ment with the corresponding analytical result, A c = 0.026, which can be ob- 
tained from(23). Only when A is increased above A c is a significant prevalence 
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observed. In Fig. 2, we plot the fraction of infected individuals in the station- 
ary state, p, for different values of /3 on scale-free network with (k) = 6. We 
observe that A c is absent. In contract with the standard SIS model, of which 
the prevalence, p, increases slowly when increasing A [24], our current epidemic 
model exhibits that p increases rapidly with A. 



6 Summary 

In this work, we analyze the SIS model that incorporates social household. We 
have focused on the impaction of geometrical property of complex networks 
and on the role of several parameters in the spreading threshold. Results show 
that the large household size N and the high within household infection rate 
are more likely to cause the spread of disease. But it's worth noticing that, even 
when local recovery rate is greater than effective infection rate, in divergent 
networks such as scale-free network, disease still can spread! This results tell 
us that even the local recover condition is good enough to give local protection, 
there are still some probability for a wide range disease spreading. It seems 
that this phenomenon can only exist in divergent networks with household 
structure. Maybe this imply that we have to care about the network structure 
much more than recover condition during disease spreading. 

Of course, the model we have studied seems more ideal. For example, we have 
supposed that the existence of the N-member households do not affect the 
property of the complex networks, and also we do not take the move of the 
individuals into account. However, the result tells us that the properties of the 
complex networks play the most important role in the epidemic spreading. 
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Captions of Figures 



Figl The fraction of the infected individuals, p,as a function of the spreading 
rate A for household structure SIS model on regular networks with k c = 4, 
N — 4. The simulations have been averaged over 200 different realizations. 

Fig2 The fraction of the infected individuals, p,as a function of the spreading A 
for household structure SIS model on scale-free networks with (k) = 6, N = 4. 
The simulations have been run in networks with 10 5 nodes. 
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Fig. 1. 



Fig. 2. 
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